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ABSTRACT 

in 

We apply our three dimensional photoionization code to model Wisconsin 
Ha Mapper observations of the H II region surrounding the 09.5V star ( Oph. 
O ' Our models investigate the porosity of the interstellar medium around £ Oph 

m , 

[N n] A6583/Ha line ratio. The ( Oph H II region has a well characterized ionizing 
^ ' source, so it is an excellent starting point for 3D models of diffuse ionized gas. We 

investigate various hierarchically clumped density structures, varying the overall 
smoothness within the clumping algorithm. By simulating the observations, we 
can estimate the porosity of the medium in the vicinity of ( Oph and find that 
within the context of our hierarchically clumped models, around 50% to 80% of 
the volume is occupied by clumps surrounded by a low density smooth medium. 
We also conclude that in order for O stars to ionize the diffuse Warm Ionized 
Medium, the O star environment must be more porous than that surrounding 
$_i ■ C Oph, with clumps occupying less than one half of the interstellar volume. 



Our clumpy models have irregular boundaries, similar to observed H II regions. 
However, in observed H II regions it is difficult to identify the precise location 
of the boundary because of the foreground and/or background emission from 
the widespread Warm Ionized Medium. This complicates the interpretation of 
the predicted rapid rise of some emission line ratios near the edge of uniform 
density H II regions and combined with the three dimensional clumpy nature of 
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the interstellar medium may explain the apparent lack of distinctive emission line 
ratios near H I H II interfaces. 

Subject headings: radiative transfer — ISM — H II regions — stars: ( Oph 

1. Introduction 

The interstellar medium (ISM) is observed to be clumpy on a very wide range of size 
scales from parsec scale molecular clouds to AU sizes inferred from interstellar scintillation 
measurements (Hill et al. 2004; Scalo & Elmegreen 2005). This dumpiness is observed to 
be self-similar over a range of scales with hierarchical or fractal dimensions in the range 
2.2 to 2.7 (Sanchez, Alfaro, & Perez 2005). The complex structures are thought to result 
from mechanical and radiative energy input into the ISM from many processes including 
stellar winds, supernovae, shocks, photoionization, and cosmic rays. This is reflected in the 
multi-phase nature of the ISM, where the hot ionized, warm ionized, warm neutral, cold 
neutral, and molecular gas co-exist. There is considerable debate among theorists regarding 
the volume filling factors of each phase and the overall structure of the ISM. Over the 
last decade large multiwavelength surveys of the ISM are providing new high spatial and 
spectral resolution observations against which observational signatures of theoretical models 
may be tested. These velocity resolved surveys include neutral gas (Hartmann & Burton 
1997), molecular gas (Dame et al. 2001), and diffuse ionized gas (Haffner et al. 2003). 
The combination of the new datasets and advances in 3D dynamical and radiation transfer 
simulations now allows for critical testing of global models of the structure and dynamics 
of the ISM. In this paper we use three-dimensional photoionization models to study the 
structure of a low density H II region in the Galaxy on scales of tens of parsecs. This is a 
first step towards understanding the structure, dumpiness, and ionization of the large scale 
diffuse ionized gas in the ISM. 

Widespread diffuse ionized gas is a significant component of the ISM in the Milky Way 
(Kulkarni & Heiles 1988; Reynolds 1995) and other galaxies (Rand 1997; Hoopes & Walterbos 
2003; Thilker et al. 2002; Collins et al. 2000). The diffuse ionized gas, hereafter referred 
to as the warm ionized medium (WIM), in the Milky Way is revealed primarily through 
faint diffuse Ha emission and pulsar dispersion measures. The physical characteristics of 
this gas, including its temperature structure and ionization state, as well as the spectrum 
of its ionizing sources are probed through analysis of optical forbidden emission lines. The 
Ha data combined with the pulsar data suggest that the WIM has a volume filling factor of 
about 20% (Reynolds 1991), and is one of the principal components of the ISM. However, 
there is considerable uncertainty regarding the structure, dynamics, source of ionization, and 
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heating of the WIM (e.g., see recent review by Reynolds, Haffher, & Madsen 2002). The 
Wisconsin Ha Mapper (WHAM) has completed an Ha survey of the Northern sky (Haffher 
et al. 2003) revealing the distribution and kinematics of the WIM in the Galaxy for the first 
time. Several large regions are also being mapped with WHAM in other lines such as H/3, 
He I A5876, [N il] A6583, [S il] A6716, [O i] A6300, and [O ill] A5007. The data have revealed 
the clumpy and filamentary structure of the WIM on large scales, such as the bipolar loop 
structure in the Perseus Arm which extends ±1 kpc above and below the Galactic plane 
(Haffner at al. 1999; Madsen 2004). In addition to the large scale structures, WHAM has 
mapped several nearby low density H II regions with sizes of tens of parsecs (e.g., Haffner et 
al. 1999). Three-dimensional photoionization modeling of these H II regions may allow us to 
estimate the structure and porosity of the ISM on these scales. Such models of well-defined 
H II regions with known ionizing sources are a first step in modeling the WIM on larger 
scales with multiple and unidentified ionizing sources. 

In this paper we study the three-dimensional structure of the H II region surrounding the 
09.5V star ( Oph, incorporating recent WHAM observations in the lines of Ha, [N il] A6583, 
and [S il] A6716 (Haffner et al. 1999; Baker et al. 2005). We use the three-dimensional 
photoionization code of Wood, Mathis, & Ercolano (2004) to model the observed Ha intensity 
and line ratios in order to estimate the dumpiness of the ISM around ( Oph. Our models 
enable us to investigate the penetration and escape of ionizing radiation from clumpy H II 
regions into the more widespread WIM. In addition, we study the observational signatures of 
ionized/neutral interfaces in smooth and clumpy models. Section 2 presents the Ha intensity 
and line ratio maps of the ( Oph region, §3 summarizes model predictions for line ratios in 
H II regions and at ionized/neutral interfaces, §4 and §5 present smooth and clumpy models 
for the observations, and we summarize our findings in §6. 



2. The C Oph H II Region 

Figure 1 shows the structure of the H II region around ( Oph. The upper left panel shows 
the Ha intensity from the SHASSA survey (Gaustad et al. 2001; Finkbeiner 2003) with an 
angular resolution of 0.8'. The upper right panel shows the same region at 1° resolution from 
the WHAM Northern Sky Survey (Haffner et al. 2003). We have interpolated the irregularly- 
spaced WHAM observations onto a regular 0?5 x 0?5 grid. The location of ( Oph is shown 
with a cross. The WHAM image shows that the overall structure of the Ha-emitting gas is 
remarkably circularly symmetric, although there are strong fluctuations on smaller angular 
scales as shown in the SHASSA image. Diffuse Ha emission associated with the ( Oph H II 
region is traced out to around 6° from the star, corresponding to a radius of about 15 pc, 
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before trailing off into more diffuse foreground/background emission toward higher Galactic 
longitudes (on the left). On the other side of the H II region the Ha intensity becomes 
confused with that from the 5 Sco H II region (see Haffner et al. 2003). We do not consider 
data from this region in our analysis. 

The sightline towards ( Oph is well studied with many absorption line determinations 
of interstellar abundances (e.g., Morton 1975; Cardelli et al. 1993, 1994; Howk & Savage 
1999). The star itself is catalogued as 09.5V, with an estimated temperature in the range 
32 000 < < 33 000 K (Code et al. 1976; Markova et al. 2004), and a distance of 140 pc 
(Perryman et al. 1997; Howk & Savage 1999). ( Oph has a heliocentric radial velocity of 
— 10.7 km s -1 , which is +3.3 km s _1 relative to the local standard of rest. The WHAM 
data shows that the ionized gas in the H II region lies at km s -1 . Its proper motion from 
HIPPARCOS is 13 ± 1 mas yr _1 in right ascension and 25.5 ± 0.5 mas yr" 1 in declination, 
resulting in a space velocity of 19.3 ±0.4 km s -1 relative to the LSR for a distance of 140 pc. 
With such a velocity, ( Oph travels more than 19 pc in 10 6 years, larger than the 15 pc 
radius of its own H II region. Therefore over its lifetime, it appears that ( Oph has moved 
from its birthplace and is now ionizing a different region of the ISM from where it was born. 
We will return to this point when discussing our photoionization models below. 

The Ha maps have not been corrected for foreground extinction by interstellar dust 
that may lie between us and the H II region. In particular, the dust lane cutting across the 
lower left quadrant in the SHASSA image is likely in the foreground and not and integral 
part of the H II region. Future H/3 observations with WHAM will allow for dust correction 
and the determination of which features are internal to the H II region and which are due to 
foreground dust. However, the [N ll]/Ha and [S ll]/Ha line ratio maps will not be affected 
by dust since Ha, [N il], and [S n] are very close in wavelength. 

In addition to the Ha maps, Fig. 1 shows WHAM line ratio maps of [N ll]/Ha and 
[S ll]/Ha. These maps show these line ratios increase with distance away from ( Oph. 
Figure 2 shows a scatter plot of Ha, [N ll]/Ha and [S ll]/Ha against radial distance from 
( Oph. The line ratios increase fairly smoothly with increasing distance and do not show 
the rapid increase expected towards the edge of a uniform density H II region (see §3). For 
comparison with models, in the rest of the paper when showing the WHAM data we show 
the mean radial surface brightness and line ratios and the dispersions about the mean. The 
mean radial surface brightness is the average value for the Ha intensity (or line ratios) in 
circular annuli centered on ( Oph. 
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3. Interface Emission: Predictions and Observations 

One dimensional photoionization models predict rapid increases in the ratios of the 
projected intensities of [O i], [N n], and [S n] relative to Ha towards the edge of uniform 
density H II regions (e.g., Henney et al. 2005). The increased line ratios occur in the 
transition region at the edge of the Stromgren sphere where the fraction of neutral gas 
is increasing and the temperature is rising due to hardening of the radiation field. The 
wavelength dependence of H° opacity allows the highest energy photons to penetrate the 
largest distances, resulting in a rising temperature in the interface region (e.g., Osterbrock 
1989). 

Observations of the diffuse ISM near the Galactic plane show low [O l]/Hct line ratios 
(Reynolds et al. 1998), suggesting that interfaces between ionized and neutral gas have 
either not been detected or do not contribute significantly to the emission. Photoionization 
models can explain the low [O iJ/Ho; observations if the WIM is highly ionized (H°/H < 
0.1 throughout), implying that ionizing radiation escapes from fully ionized (i.e., density 
bounded) H II regions having few H I H II interfaces (e.g., Mathis 1986; Domgoergen & 
Mathis 1994; Mathis 2000; Sembach et al. 2000). Also, recent work by Giammanco et al. 
(2004) showed that if a clumpy H II region is comprised of a mixture of fully and partially 
ionized spherical clouds, the line ratios increase more slowly with radius than in uniform 
density, ionization bounded models (e.g., see the radial dependence of [O i]/H/3 in their 
Fig. 6). This is because the combination of many spherical clouds along a given sightline 
through the H II region sample a wide range of ionization parameters and hence the line 
ratios do not show the rapid increase of uniform, ionization bounded models. 

While the predicted rapid rise of temperatures and line ratios at the edges of H II regions 
has not been observed (e.g., Pauls & Wilson 1977), gradual rises in the line ratios have been 
observed with height above the midplane ionizing sources in the Milky Way (Haffner et 
al. 1999) and several other edge-on galaxies (Rand 1998; Otte et al. 2001, 2002). The 
elevated [N nj/Hoj line ratios observed at large distances from the midplane in the Perseus 
Arm (Haffner et al. 1999) have been interpreted as due to an extra heat source in addition 
to photoionization heating (Reynolds et al. 1999). This additional heat source must have 
the property that it dominates over photoionization heating at low densities. A constant 
energy source or one that is proportional to n c would suffice, since photoionization heating 
is proportional to n 2 e . Candidates for the additional heating include photoelectric heating 
from grains (Reynolds & Cox 1992), magnetic reconnection (Raymond 1992), dissipation of 
turbulence (Minter & Spangler 1997), shocks and cooling hot gas (Slavin, Shull, & Begelman 
1993; Collins & Rand 2001). However, hardening of the radiation field also increases the 
temperature of photoionized gas and may help to explain in part some of the observed 
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elevated line ratios (e.g., Bland-Hawthorn, Freeman, & Quinn 1997; Elwert 2003; Wood & 
Mathis 2004). 

In this paper we present models for the H II region associated with ( Oph. Additional 
heating will not be important for this source since the gas density is larger and the line ratios 
are much lower than observed at high latitude in the Perseus arm. The WHAM observations 
presented in Figures 1 and 2 show that [N ll]/Hct increases away from the source, but does 
not appear to show the rapid rise at the edge of the H II region predicted by spherical models 
with uniform density (see §4.1 below). To understand this behavior we have investigated 
photoionization models that incorporate various smooth and clumpy density structures for 
the H II region and the effects of diffuse foreground/background emission. We discuss the 
plausibility of the smooth models and the insight three-dimensional models can provide into 
the structure of the ISM around ( Oph and the escape of ionizing photons from H II regions. 
Our 3D models produce similar results to those of Giammanco et al. (2004) for radial line 
ratio gradients, but our 3D radiation transfer allows for shadowing of clumps and multiple 
clumps along any sightline from the star. 



4. Photoionization Models 

To model the WHAM observations of the ( Oph H II region, we use the three-dimensional 
photoionization code described in Wood, Mathis, & Ercolano (2004). This code calculates 
the 3D ionization and temperature structure for arbitrary geometries and illuminations, 
keeping track of the ionization structure of H, He, C, N, O, Ne, and S. The opacity is from 
to H° and He and we only consider photon energies in the range 13.6 eV to 54 eV. Heating 
is from photoionization of H and He and cooling is from H and He recombination, free-free 
emission, and collisionally excited line emission from C, N, O, Ne, and S. The output of our 
code is the ionization and temperature structure from which we calculate emissivities and 
intensity maps for the various emission lines we wish to study. A complete description of 
the code is presented in Wood et al. (2004). 

We perform the radiation transfer for stellar and diffuse recombination ionizing radiation 
on a 65 3 linear Cartesian grid. Due to the resolution of our grid, we do not consider emission 
from cells that are more neutral than H°/H > 0.25 (see discussion in Wood et al. 2004). 
These cells lie towards the edge of the ionized volume and account for less than 20% of the 
ionized cells in our ionization bounded models. We have constructed intensity and line ratio 
maps including these cells and find that our results for the [N ii]/Hq; line ratio maps do 
not change at the ~ 5% level. Moreover, in the cells at the ionization boundary our code 
is incomplete, since we ignore the effects of dynamics and shocks at the ionization front 
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(Henney et al. 2005). Higher resolution simulations (128 3 grids) do not appreciably change 
our results, again at around the ~ 5% level. 

Our code performs well compared to other photoionization codes in predicting ionization 
fractions, temperatures, and line strengths. Here we focus on modeling the observed Ha 
intensity and [N ll]/Ha maps and do not consider the [S ll]/Ha observations. Modeling 
emission from sulfur is problematic due to the unknown dielectronic recombination rates for 
S (Ali et al. 1990). In addition we do not consider the effects of dust within the H II region 
and note that our temperatures may be slightly higher than calculations that include cooling 
from more elements (e.g., see Sembach et al. 2000). Since the goal is to investigate the 3D 
structure of the H II region by modeling variations in the line ratios towards the edge of the 
H II region, ignoring these effects will not change our conclusions on the H II region density 
structure. Inclusion of these additional effects will be required for more comprehensive 
models of future observations of H/3, [O i], [O il], [O ill], and He I. In particular, we may 
be able to model data from several H II regions to place constraints on the S dielectronic 
recombination rates. 

Relative to H, the adopted abundances for He, C, N, O, Ne, and S in our models 
are 0.1, 140 ppm, 75 ppm, 319 ppm, 117 ppm, and 18 ppm. These are interstellar gas 
phase abundances from the compilation of Sembach et al. (2000). The ionizing spectrum is 
taken to be that of a 32 000 K WM-basic model atmosphere from the library of Sternberg, 
Hoffmann, & Pauldrach (2003). This temperature is within the range estimated for an 
09.5V star (see compilation of temperature/spectral type scales and discussion in Harries, 
Hilditch, & Howarth 2003, Table 4). The WM-basic model atmospheres include the effects 
of line blanketing and stellar winds (Pauldrach et al. 2001). The ionizing luminosity and 
density are varied in the models to reproduce the observations of Ha and [N ll]/Ha. The 
following sections present results for smooth models (with and without density gradients) 
and three-dimensional, hierarchically clumped models. 



Our spherically symmetric models are centered on a source with an ionizing photon 
luminosity Q (s _1 ), and the density within the cloud (cm" 3 ) is given by 



where the parameter p controls the radial density gradient within the cloud. In all models 
(smooth and clumpy) we have left the inner 10% of the grid empty, which prevents the 



4.1. Smooth Spherically Symmetric Models 




(1) 
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clumping algorithm randomly placing the star in the center of a dense clump. The grid is 
60 pc on a side and contains 65 3 cells, while the observed ( Oph H II region is about 40 pc in 
diameter but has some extensions to greater distances. We have investigated many (p, n , 
Q) combinations for ionization and density bounded models and show those that reproduce 
the overall Ha intensity level observed by WHAM. Table 1 summarizes the range of adopted 
model parameters and the results of the models are compared to the observations in figure 3. 

The radial variation of Ha and [N ll]/Ha may be influenced by the presence of fore- 
ground and background emission. This is illustrated in Figure 3 which shows models without 
(left panels) and with (right panels) the addition of a constant intensity representing fore- 
ground/background emission. Both the smooth (shown here) and clumpy models (shown in 
§4.2) have been convolved with a Gaussian beam to simulate the 1° angular resolution of 
WHAM. The diffuse background is taken to be 1.5 R at Ha and 0.6 R for the [N il] simu- 
lation. These values provide a good fit to the observations (see Fig. 2) and also are typical 
for the diffuse interstellar medium near the Galactic latitude (b = 24°) of ( Oph (see also 
Haffner et al. 2003; Madsen 2004). The addition of diffuse foreground/background emission 
(which is an unavoidable contaminant in real observations) clearly suppresses the very large 
line ratios that are predicted in the faintest outermost regions (~ 20 pc) of the H II region. 
However, even with the addition of foreground/background emission, the uniform density, 
ionization bounded model (i.e., a classic Stromgren sphere) exhibits a rise in the [N ll]/Ha at 
the edge of the H II region that is steeper than the observations. The Ha intensity gradient 
in such models is also much steeper at the edge of the ionized region than observed (see 
also the giant H II region models of Giammanco et al. 2004). Density bounded or "leaky 
H II region" models (Fig. 3) do not exhibit the rapid rise in line ratios because there is no 
ionized/neutral interface in the simulations (e.g., Mathis 2000; Sembach et al. 2000). The 
density bounded models shown in Fig. 3 allow about 40% of the ionizing photons to escape. 
Changing the radial density gradient from uniform (p = 0) to steeper laws (p — 1,2) changes 
the Ha radial surface brightness profile. Models with gradients steeper than p = 1 do not 
match the WHAM data. Models with p = 1 appear to give a reasonable match to the Ha 
intensity, but all the density bounded models underpredict the [N ll]/Ha observations. 

As mentioned previously, ( Oph is at a different radial velocity from its associated H II 
region, so models that place the star at the centre of a spherical cloud with radial density 
gradients seem contrived. A more realistic scenario is that ( Oph is ionizing a clumpy 
medium and we explore 3D cloud geometries and the effects on the intensity and line ratio 
maps in the next section. 
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4.2. Hierarchically Clumped Models 

There is much evidence for hierarchical structure in the ISM, with surveys of clouds 
revealing fractal structures (e.g., Elmegreen & Falgarone 1996). We therefore investigate 
3D H II regions which have hierarchically clumped density structures. A 3D hierarchical 
density structure is created using the algorithm presented by Elmegreen (1997). We use 
five hierarchical levels and the density in the grid is set as follows. At the first level we 
randomly cast iV points with x, y, and z coordinates in the range (0,1). At each subsequent 
hierarchical level we cast iV random points around each of the points cast at the previous 
level. The casting length for the x, y, and z points at each subsequent level is in the range 
±/\( 1 ~ H ) /2, where H is the hierarchical level being cast, and / = log N/ log A is the fractal 
dimension of the hierarchical structure. The casting length, A, gets smaller for subsequent 
hierarchical levels. Note that in this description A is dimensionless and we construct the 
hierarchical grid within a cube with side of unit length. Actual dimensions are constructed 
by mapping the hierarchical grid onto our 3D density grid, so that the unit cube corresponds 
to the physical size of our density grid. We use a fractal dimension / = 2.6, appropriate for 
interstellar clouds (e.g., Sanchez, Alfaro, & Perez 2005). Using this algorithm, the density of 
a cell in our grid is proportional to the number of points cast at the final level that lie in the 
cell. If points are cast beyond the x, y, or z boundaries of the grid, they are added into the 
corresponding cells on the opposite side of the grid, e.g. if a point is cast at x — x max + e, we 
change its coordinate to be x = e and similarly for y and z coordinates. Figure 4 illustrates 
this algorithm in 2D for a three tier hierarchical scheme, with four random points cast at 
each level. The filled in squares are the first castings, diamonds the second level, and crosses 
at the third level. The dotted lines show overlaid grid cells. We assume that a fraction, 
/smooth, of the average density is present in every cell. The remaining mass is distributed in 
proportion to the number of crosses that are in each cell. 

The smooth component represents unresolved fractal structure that cannot be resolved 
within our grid. In real nebulae the hot, shocked stellar winds may create very low density 
regions between the clumps of gas. The resulting clumpy density grid is renormalized so that 
it has the same mean density structure as the smooth grid. We generally adopt /smooth = 1/3, 
but also consider the case where there is zero smooth density component. Our adopted value 
of /smooth comes from radiation transfer models of the penetration of radiation into turbulent 
cloud models from hydrodynamic simulations. A value of / sm ooth — 1/3 provides a good 
match for the internal intensity levels between the hydrodynamical simulations and density 
structures generated with the fractal algorithm (Bethell et al. 2004). 

To investigate the dumpiness of the ISM around ( Oph we varied iV at the first hierar- 
chical level, using Aq = 32, 64, 128, 256, and 512, but keeping N = 32 at all four subsequent 
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hierarchical levels. With this modification to the clumping algorithm the relation between 
casting length, A, and fractal dimension, / = log N/ log A, is determined assuming N = 32. 
Increasing N\ at the first level gives a progressively smoother density structure at large spa- 
tial scales, but maintains the observed fractal dimension of interstellar clouds on smaller 
scales. We have not investigated the effects on the ionization structure and line ratios of 
changing the value N = 32 at the second and higher hierarchical levels. Increasing N will 
give a larger dynamic range and finer structure in the hierarchical density grid, because the 
total number of points used to determine the density is given by the number of castings 
at the final level, N\ x N 4 . The choice of N = 32 is not observationally motivated, but is 
determined by the resolution of our density grid. Using very large values of N is not justified 
because the fine structure will not be resolved by our grid. 

The algorithm produces maximum density contrasts between the densest clumps and 
the smooth component of 65, 35, 20, 10, for N\ = 32, 64, 128, and 512 respectively. A small 
Ni allows ionizing radiation to penetrate further between the clumps than in a uniform 
medium (Elmegreen 1997). The fraction of cells that contain only the smooth component of 
density, without any of the points that were cast randomly, is 65%, 45%, 20%, and 0.1% for 
these models. If the density in the smooth component is almost zero, ionizing radiation can 
escape along the empty sightlines into the larger scale ISM. 

Having a smooth component of zero density did not significantly change the morphology 
of our Ha maps compared to the /smooth — 1/3 that we adopt because the emission is 
dominated by the higher density ionized clumps. For iV\ = 32 and / smoo th = 1/3, the smooth 
component accounts for about 1/3 of the Ha intensity level. The smooth contribution to the 
Ha intensity decreases with increasing N±. Elmegreen (1997) speculates that much of the 
volume of a fractal ISM is at very low density, like our models with no smooth component. 

We investigated many different clumpy models for each Ni value by changing the random 
number seed used for setting up the hierarchical density grid. For a given N\, each model has 
a different density structure and hence slightly different Ha and [N il] intensity profile due to 
the random casting of clumps around the star, but the overall features in the intensity maps 
are similar. For each Ni of our clumpy models, Figures 5 through 8 show one pixel thick slices 
through one of the density grids and slices showing the hydrogen ionization fraction H°/H. 
The figures also show the Ha intensity and [N ll]/Ha line ratio maps, without the addition 
of uniform foreground/background emission. Figures 9 to 12 show the radial variation of Ha 
and [N ll]/Ha with the addition of uniform foreground/background emission. For each Ni 
model, the figures show one model with the standard deviations about the mean and also 
the radial profiles for simulations of five different density grids. 
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5. Results 

Our models are constrained by the emission line intensities observed at various angles 
from the star (or distances projected onto the sky) and their dispersions. For all clumpy 
models the ionizing luminosity is Q = 8 x 10 47 s _1 and the mean density of the H II region 
is n = 2cm~ 3 . Both the luminosity and density are about 50% smaller than estimated by 
Elmergreen (1975) from analysis of older Ha data from Reynolds et al. (1974). Our estimates 
for the luminosity and density are well constrained by the intensity and angular extent of 
the Ha data. We find that N\ ~ 64, followed by coarser clumping on the smaller scales, 
provides the best fit to the WHAM data. This implies that the real ISM in this direction is 
not hierarchical, but is smoother on a scale of ~ 10 pc than at smaller scales. This result 
reflects the general circular symmetry in the WHAM Ha image in Fig. 1, but the small scale 
irregularities are not revealed in the WHAM data because of its relatively coarse angular 
resolution (1°). 

As expected, Figures 5 through 12 show that the clumpy model that most resembles 
a Stromgren sphere has Ni = 512, since distributing many random points approximates 
a uniform distribution. This model has a somewhat irregular boundary, but there are no 
sightlines from the star that do not intersect dense clumps, so overall the model has a very 
smooth structure and exhibits the large increase in the [N ll]/Ha line ratio at the edge of the 
H II region. As Ni is reduced, the medium becomes less smooth and there are sightlines from 
the star through which ionizing radiation can traverse large distances before intersecting a 
clump, or not intersect any clumps. Models with N\ = 32 show a gradual radial decline of 
the Ha intensity. As N 1 is increased, the fraction of low density sightlines decreases and 
the Ha exhibits a steeper radial gradient. The concavity of the radial Ha surface brightness 
profiles change for larger N 1 and approach that of the uniform density Stromgren sphere 
in Figure 3. The trends in our models are similar to the radial surface brightness profiles 
and line ratios found by Giammanco et al. (2004) in their H II region models comprising 
spherical blobs surrounding an ionizing source. 

Within the context of our clumping algorithm, we find that models with Ni = 32 
produce too much Ha at large radii and variations about the mean radial intensity are 
larger than observed. Models with Ni > 256 resemble Stromgren spheres, producing too 
steep a gradient of the Ha and a rapid rise of [N Il]/Ha at the edge of the H II region 
and small variations about the mean values. Models with Ni = 64 and 128 provide better 
matches to the WHAM data in reproducing the observed levels of Ha and [N ll]/Ha, their 
radial variations, and standard deviations. 

While the intensity maps for this low density H II region model are not very sensitive to 
the smooth component, the smooth component has a large influence on the ionizing radiation 
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that may escape the H II region and ionize the larger scale diffuse ionized gas. Models with 
Ni = 32, 64, and 128 and zero smooth density component (/smooth — 0) typically allow around 
10%-25%, 3%-15%, and l%-3% respectively of the stellar ionizing photons to escape the 
H II region directly along the zero density paths from the star. The large variations in the 
escape fractions for each iVi value are due to the random placing of clumps around the star. 
Our 3D photoionization models therefore confirm the analysis of Elmegreen (1997) that a 
fractal ISM will allow for ionizing photons to traverse much larger distances than a smooth 
density structure. Further clearing of a fractal ISM by feedback from photoionization (Dale 
et al. 2005) and stellar winds will likely provide an even larger escape fraction from H II 
regions than those determined from our models. 

Since the widespread Ha emission from the Galactic WIM requires about 15% of the 
ionizing photons from Galactic O stars (Reynolds 1995), it appears that a 3D fractal ISM 
structure could allow for such leakage from H II regions. Our hierarchical density models 
for the ISM around ( Oph with Ni = 64 and 128 typically provide escape fractions of less 
than 15%, and even lower if the density in the smooth component is not zero. Therefore, the 
ionizing sources responsible for powering the WIM in the Galaxy must reside in regions of 
the Galaxy that are more porous than the region around £ Oph. Perhaps their environment 
has been shaped more by turbulent motions or are regions cleared out by stellar winds, thus 
enabling the escape of ionizing photons. Future models of WHAM observations of other H II 
regions will enable us to determine whether the porosity levels for the ISM around ( Oph 
are typical. 

Our models placed the star at the center of a fractal cloud, so models where the star 
is at the edge of the cloud will allow for larger leakage, since the cloud no longer covers the 
whole sky as seen from the star. The major question from this work relating to ionization of 
the diffuse ionized gas is what is the density of the smooth component in fractal models? We 
cannot determine how much Ha emission from the WIM comes from an extensive, relatively 
smooth component and how much is from ionized clumps and the ionized faces of clouds. 
Further progress to answer this question may be made by conducting photoionization simu- 
lations on the 3D density grids being produced by global models of the ISM (e.g., de Avillez 
& Berry 2001). Such models predict the ISM density structure and the locations of ionizing 
sources and their validity may be tested against observations through 3D photoionization 
simulations. 
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6. Summary 

We have presented smooth and clumpy 3D photoionization models for the ( Oph H II 
region. We find that a star with T = 32 000 K and ionizing photon luminosity Q = 8 x 10 47 s" 1 
and mean density of the H II region of n ~ 2 cm -3 reproduce the observations. These numbers 
are well constrained by the WHAM data. The stellar temperature and ionizing luminosity 
are remarkably close to recent determinations of stellar properties for 09.5V stars (Martins 
et al. 2005, Table 4). Larger densities require a larger ionizing luminosity to match the 
extent of the H II region, but then the Ha values are larger than observed. Smaller densities 
have the opposite effect, requiring a lower ionizing luminosity to match the extent of the H II 
region, but then the simulated Ha intensity is much smaller than observed. 

Smooth models that best match the observations are density bounded and have a radial 
density gradient that falls off no faster than r~ l . However, the velocity offset between ( Oph 
and its H II region suggest that these uniform models are unrealistic as they require the chance 
alignment of the star at the center of the cloud. The inclusion of foreground/background 
emission suppresses the rapid rise in the [N ll]/Ha line ratio predicted at the edge of uniform 
Stromgren spheres. 

Hierarchically clumped models, which reproduce the observed fractal structure in the 
ISM, also reproduce the Ha surface brightness and the shallow gradient of [N ll]/Ha away 
from ( Oph. The best models suggest that around 20% (Aq = 128) to 50% (Aq = 64) of 
the H II region volume is occupied by gas distributed in a smooth component or at very low 
density, with the remainder in clumps. Three-dimensional models that are more uniform do 
not match the observations. 

The large, low density regions in 3D models could allow ionizing photons to escape 
at levels that sustain the ionization of the large scale diffuse ionized gas in the Galaxy. 
However, the escape fractions from our preferred clumpy models for ( Oph are less than 
the 15% escape fraction required to ionize the WIM. Therefore we conclude that the stars 
responsible for ionizing the WIM must reside in regions of the ISM that are more porous than 
that surrounding ( Oph, possibly evacuated by supernovae and/or stellar winds. A more 
detailed investigation into the leakage of ionizing photons to the WIM requires consideration 
of dynamical clearing and the time evolution of the ionizing luminosity of O stars. 

Future work will include more comprehensive modeling of the planned mapping obser- 
vations of ( Oph in additional emission lines. The WHAM survey has mapped out several 
other low emission measure H II regions and 3D modeling will allow us to test whether the 
porosity values we estimate for the ( Oph region are common in other regions. Modeling the 
[S II] emission from many such regions may help to place constraints on the as yet unknown 
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dielectronic recombination rates for sulfur. From a theoretical perspective, global dynami- 
cal models of the ISM can be tested by applying 3D photoionization codes and comparing 
their observational signatures with the high resolution observations now available of diffuse 
ionized gas in our own and other galaxies. 

We thank Lynn Matthews and Jane Greaves for comments on an early version of this 
paper. We acknowledge funding from a PPARC Advanced Fellowship (KW) and a PPARC 
Visiting Fellowship to the University of St Andrews (JSM, RJR, GM). WHAM is funded by 
grants from the NSF (AST 0204973) and the University of Wisconsin Graduate School. 
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Table 1: Smooth Model Parameters 



n 
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^max 
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Ionization 


(cm -3 ) 
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bounded? 
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no 
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Fig. 1. — Views of the ionized gas near ( Oph (Sharpless 27). The top panels display the loga- 
rithm of the Ha emission from (left) a moderate-resolution (~ 8') image provided by SHASSA 
(Gaustad et al. 2001; Finkbeiner 2003) and (right) a lower-resolution representation from 
WHAM. The WHAM spectra have been integrated between —25 km s _1 < t>LSR < +25 
km s _1 for these images. The lower panels show the (left) [N II]/Ha and (right) [S II]/Ha 
ratios of the region. The location of £ Oph is denoted by the '+' in each panel. Part of 
the S Sco H II region may be seen at the right edge of the SHASSA image. All four images 
are displayed with an Aitoff-Hammer projection. The axis labels are Galactic co-ordinates 
(latitude vs. longitude) in units of degrees. 
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Fig. 2. — WHAM data showing the variation of Ha, [N n]/Ha, and [S n]/Ha as a function 
of radius from ( Oph. Notice the relatively smooth increase of the line ratios with distance 
and the absence of the rapid increase expected at the edge of a uniform density Stromgren 
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Fig. 3. — Smooth models without (left) and with (right) the addition of fore- 
ground/background emission. The heavy solid line shows the mean radial intensity and 
standard deviation (error bars) measured by WHAM. The lighter lines show the models 
from Table 1: uniform ionization bounded (solid), density bounded: uniform (dotted), 1/r 
density (dashed), and 1/r 2 density (dot-dashed). 



-22 - 



.0 



0.8 



0.6 



0.4 



0.2 



0.0 
0.0 



1 1 

— 


1 1 

+ 


*■ + 




- 


+ 


-+-+ 

/\ 1 

° •> 

o + 

+ + 


+ 
+ 


--4 


— 


+ 
+ 


.V: 




+ + + + 

o ■ 

+ 

++ 


i - 


o 
+c 

+ + 


+ 
+ 
>+ 


+ 


- - -o- ----- 

+ 

+ 

■ + + 




1 1 


i i 


+ 


rr#<?< 

+ 

+ + + 


> 

+ 



0.2 



0.4 0.6 

x 



0. 



1 .0 



Fig. 4. — A 2D illustration of the fractal generating algorithm for a three tier scheme. Four 
random points (filled in squares) are cast at the first level. At each subsequent level a further 
four points are cast around each of the points cast at the previous level. In this example there 
are a total of 16 points (diamonds) cast at the second level, and a total of 64 points (crosses) 
cast at the third level. The density in each cell in the grid (dashed lines) is proportional to 
the number of crosses in the cell. In our 3D models we use a five tier scheme with Aq = 32, 
64, 128, or 512 points cast at the first level and at each subsequent level we cast N = 32 
points around each of the points cast at the previous level. Increasing Aq results in a overall 
smoother medium (see Figs. 5 through 8). 
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Fig. 5. — Hierarchically clumped model with N\ = 32 (see text) showing slices through the 
center of our grid of the number density and hydrogen ionization fraction (upper panels) and 
the projected Ha intensity and projected [N ll]/Ha line ratio maps (lower panels). Axes are 
labeled in parsecs. This model allows for many sightlines from the star that do not intersect 
a dense clump, producing extended low intensity Ha emission which may get lost in strong 
foreground/background emission. 
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Fig. 6.— As Fig. 4, but for JVi = 64. 
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Fig. 7.— As Fig. 4, but for JVi = 128. 
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Fig. 8. — As Fig. 4, but for Ni = 512. Notice how much more circularly symmetric the 
ionized region (upper right) is than in Fig. 5, because the many points in the initial casting 
tend to fill space unifromly. 
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Fig. 9. — The Ha intensity and [N ll]/Ha line ratios for the clumpy models with N% = 32, 
plotted against the projected offset distance from ( Oph. Left panels: red lines show the 
mean radial intensity and standard deviations measured by WHAM and black lines show the 
same for one of the N% = 32 models. A constant intensity of 1.5 R and 0.6 R has been added 
to the model Ha and [N n] maps to simulate the diffuse foreground/background emission 
around ( Oph. Right panels: red lines again show the WHAM data and the black lines show 
five different N± = 32 models — each model has a different randomly cast hierarchical density 
structure. These models show high Ha at large radii due to ionizing photons propagating to 
large distances through the very porous density structure (see Fig. 5). The [N ll]/Ha ratio 
is systematically lower than observed and the standard deviations of the models are larger 
than in the WHAM data. 
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Fig. 12. — As Fig. 8, but for N\ = 512. This clumpy model closely resembles the uniform 
density, ionization bounded H II region since this model has the smoothest density structure 
of the 3D models presented (see Figs. 5 through 8). 



